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ABSTRACT 


In studying the kinematics of the formation, reactions, and de- 


o and St Sh 


cay of unstable strange particles such as ae yx a tc] 
a liquid hydrogen bubble chamber a range-energy relation is essential. 
To obtain such a relation, one must know the density of superheated 
liquid parahydrogen. Heretofore, this density has not been determined. 
Direct measurement by conventional means of the density of the liquid 
hydrogen in the bubble chamber is difficult because the liquid is super- 
heated for only 20 to 30 milliseconds. 

To determine the density, the mean range of ur mesons produced 
by 1’ mesons decaying at rest in the 10-inch UCRL liquid hydrogen 
bubble chamber was measured and found to be 1.103 + 0.003 cm 
(standard deviation 5.0%; 444 1 > u-—-e events). From the Bethe-Bloch 
theory anda Ti kinetic energy of 4.12 + 0.02 Mev., this value of range 
gave a density of (5.83 + .06) x ow By oe for superheated liquid para - 
hydrogen. The temperature and pressure of the hydrogen during the 
sensitive time were 27.6 + 0.1°K and 48+ 5 psia. A plot of the range - 
energy relation for protons, tT, yw, K, 5, and == particles based on 
Bethe-Bloch theory, and the above value of density is given for kinetic 


energy values ranging from 1 to 1000 Mev. 
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l. Introduction 


Since the first successful operation of a bubble chamber by 
Glaser [1] in 1952, its use as an instrument for nuclear research has 
grown rapidly and the information yielded has aided our understanding 
of nuclear physics. At the present time various types of chambers are 
in use or under development. These differ in many respects, but 
primarily they can be classified by the type of sensitive material used, 
1.e., hydrogen, hydrocarbons, or--more recently--helium. The 
advantages and disadvantages of using hydrogen as the active material 
are fully discussed elsewhere [2]. At the University of California 
Radiation Laboratory a 10-inch hydrogen bubble chamber is being used 
in the study of the production and decay of short-lived unstable particles 
such as i, oP, and = . An accurate range-energy relation is 
needed in determining the kinematics of the processes involved. If 
the density of the liquid hydrogen were known, one could use the Bethe - 
Bloch theory [3,4] to determine the range-energy relation. It is, 
however, very difficult to measure the density of superheated liquid 
hydrogen, the condition in which it is sensitive to ionizing particles, 
because it remains in this condition only 20 to 30 msec. One could 
extrapolate into the superheated region, using measurements made on 
the normal liquid, but the accuracy of this method is unsatisfactory for 
studies of the unstable particles mentioned above. 

One point on the range-energy relation can be found if we measure 
the range of a particle of known initial and final kinetic energy. Using 
this experimental point we can calculate the density of hydrogen for use 
in the Bethe-Bloch theory to determine the range -energy relation over 
a wide range of energies. 

The method selected for the experimental determination of a point 
on the range -energy curve was to measure the range of the ut meson 


i ut ie decay chain. The selection was based on such 


in the T 
factors as the accuracy with which the ra energy is known, the 


accuracy with which the ra range could be measured, and the 
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availability of 1’ mesons in sufficient quantity ald at the preper Ghergy 
range. 

The predominant mode of decay of the a’ mesonis wv’ > ut + Vv. 
Experimental evidence for other modes of decay for the Tr’ exists Lode 
but the branching ratios are so small aioe to oe) that for our 
purposes the decay of the mw’ can be considered a two-body decay yield- 
ing a ra with 4.12 + 0.02 Mev [6] initial kinetic energy if the decay 
takes place at rest. 

The ut is known to decay at rest according to the relationship 
ut > et +a Gav. 

The e’ has a finite energy spectrum, because the decay involves more 
than two products. The important feature of the ut decay is that it 
does take place at rest, i.e. at zero kinetic energy, andis evidenced 
by the formation of the positron. Knowing the initial energy of the ur, 
we need only measure its range in the bubble chamber to obtain a point 
on the range energy curve, and as explained above, the range -energy 
relation. In the actual extrapolation of the range-energy relation from 
the point determined to other energies itis simpler to use the data on 
protons available in the literature. Aron [7] has computed range-energy 
tables for protons in hydrogen at standard conditions. These tables 

are based on a theoretically determined mean-excitation potential of 
17.5 ev [8] for hydrogen. By properly scaling the data we can calculate 
the range -energy relation for hydrogen under the operating conditions 
of the bubble chamber. 

Since the kinetic energy of the i at the time of its formation is 
known to only +1/2 % [6], a final range error of +1/3% seemed a 
reasonable objective. The factors affecting the final range error are 
range straggling and measurement errors. If measurement errors can 
be made small compared to the range straggling, the accuracy of our 
range determination will be limited by the number N of ur ranges 


measured, 1.¢€., 
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Then for a measurement plus straggling error of 5%, we need to 


range error ¥ (straggling + measurement error)/N 


measure about 300 ut ranges to reduce our range error to 1/3%. 

To minimize the measurement errors, certain restrictions in the 
selection of 7 > ur +e’ events for measurement had to be observed. 
Since tracks in the chamber consist of a series of bubbles separated by 
finite distances, one cannot assume that the first and last bubble in any 
series marks the beginning and end of the particle track. In 
1 > ra ee" decay events, the points of titers ection between the T 
and the uw and between the e andthe yp determine the limits of the 
wy. travel. These points of intersection, once determined, offer the 
easiest method of measuring the length of the w track. The error in 
the determination of a point of intersection is inversely proportional 
to the sine of the angle of intersection, and both angles of intersection 
enter into the final range accuracy. 

As may be seenin Fig. 1, the stereo angle of the camera used 
with the 10-inch bubble chamber is ~4°. This small angle makes 
measurement errors sensitive to the error in depth determination. 
Specifically, the error in range is proportional to the product of the 
error in the measurement of change in depth of the particle track, and 
the sine of the angle of dip of the particle track, i.e. measurement 
error =(errorin Az) sink. 

Finally, an error in the measurement of range may be introduced 
by 1 mesons that decay in flight but cannot be recognized as doing so. 
On the basis of a wm lifetime of 7 = 2.56 x co sec [9], the probability 
of decay in flight in the last centimeter of range is calculated to be 
0.0076. This error is obviously insignificant. The o, being much 
longer-lived [9] that the nv , has a correspondingly smaller chance 
to decay in flight. 

With the range-energy relation determined, and the density of the 
superheated liquid hydrogen calculated, we can comment on the ad- 


visability of extrapolating the known data on normal liquid hydrogen 


into the superheated region. 
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Fig. 1. Line drawing of 10-inch hydrogen bubble chamber 





2. Experimental Arrangement 


The 1’ beam was obtained by bombarding the 79° lead target of the 
Bevatron with protons of 3 Bev energy; tT’ mesons of approximately 
245 Mev/c momentum were emitted in the backward direction and 
deflected out through the 76° window by the Bevatron magnetic field. 

A small ''C'' magnet placed at the exit window was used to make minor 
corrections in the directions of the m’ beam. After being focused by a 
4-inch quadrupole magnet, the beam was deflected through 60° and 
collimated into a 4-by-4-inch cross section for entry into the 10-inch 
bubble chamber. The quadrupole -magnet current was adjusted to 

focus the beam at the bubble -chamber-entrance window. A schematic 
of the experimental arrangement is shown in Fig. 2. 

To stop the r mesons in the center of the chamber, the beam 
momentum had to be reduced from 245 to 120 Mev/c prior to entry 
of the beam into the vacuum jacket surrounding the chamber. The 
value 120 Mev/c takes into consideration the stopping power of the 
vacuum jacket, the chamber window, and 5 inches (chamber radius) 
of liquid hydrogen at a density of about 0. 06 g/cm”. The true energy 
of the 1’ beam was determined by taking a range curve with copper 
as the absorber. The Bevatron beam intensity was 2 x 107 protons per 
pulse. This was calculated to give an average of three wT" per pulse 
stopping in the chamber. Some of the pulses yielded many more, as 
shown in Fig. 3. 

Although the chamber was surrounded with lead shielding, it was 
necessary to turn on the bubble chamber magnet to reduce the back- 
ground of low-energy particles, mostly electrons, entering the chamber 
and obscuring the decay events. The curvature of charged-particle 
tracks caused by the magnetic field of the chamber introduces a 
systematic error in the range measurements. The error introduced 
was found to be insignificant (see discussion of magnetic curvature in 


Section 5E). 
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Locating, positioning, and focusing the beam was done by use of 
coincidence techniques with scintillation counters. Background effects 
were reduced by gating the scalers to count only during the Bevatron 


pulse, 
3. Measurement of Liquid Hydrogen Temperature and Pressure 
A. Temperature 


Temperature measurements in the 10-inch bubble chamber are 
made with two hydrogen-filled vapor-pressure thermometers whose bulbs 
are located one above the other about 3.6 inches apart on the inside wall 
of the chamber (see Fig. 4). The bulbs are connected to opposite sides 
of a pressure-difference gauge located outside the chamber. This gauge 
provides a meansfor measuring the vertical temperature gradient during 
the operation of the chamber. The lower of the two bulbs is also con- 
nected to a Bourdon pressure gauge. 

During chamber operation the temperature is automatically controlled 
with heaters surrounding the chamber wall [10]. The pressure in the 
lower bulb is used in the controlling. The desired chamber temperature 
is preset on the controller, and this temperature is maintained within the 
accuracy of the Bourdon Gauge (+ 0.2 psig = + 0.02 °K at 65.5 psig). The 
vapor-pressure difference varies slightly during operation, The chamber 
temperature was checked continuously and recorded at least once every 5 
to 10 minutes during the 3 hours the experiment lasted. The reading of 
the lower vapor-pressure thermometer remained essentially constant, 
and equal to the value set on the Bourdon Gauge. According to the experi- 
mentally determined vapor-pressure curve for liquid hydrogen[11], the 
65.5-psig reading of the lower thermometer is equivalent to 27.8 °K, and 
the average vapor-pressure difference of 3.5 psi corresponds to a tempera- 


ture gradient of 0.021 °K/cm. 
B. Pressure 


The problem of determining the pressure to which the hydrogen ex- 


pands was far more difficult than the temperature measurements. Built 
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Fig. 4. Sketch of v ane -pr e-—-temperature bulbs in the 10-inch 
bubble chambe 





into the walls of the chamber is a pressure-sensing device referred to 
locally as a Linlor pressure gauge [12]. The gauge gives an indication 
of the instantaneous pressure in the chamber, and is used primarily in 
adjusting the bubble chamber operating cycle. The gauge basically 
operates as follows: a change in pressure in the chamber causes a 
change in capacitance in an electrical circuit, which is converted into 6 
change in voltage and displayed on an oscilloscope. In principle the 
method is quite simple, and the gauge accomplishes the purpose for which 
it was designed. However, the gauge does have certain limitations that 
could not be eliminated. The limitations were: (a) the sensitivity of the 
gauge was very low, (b). the reference voltage of the onaee as well as 

the sensitivity was subject to drift, and (c) the response of the gauge 

was nonlinear. 

The pressure measurements were made in the following way. With 
each Bevatron pulse, the voltage output signal from the Linlor gauge 
was displayed on an oscilloscope. The oscilloscope screen was monitored 
continuously, and a series of pictures taken every 5 to 10 minutes to 
correspond with temperature measurements. Land camera pictures 
gave better results than a 35-mm scope camera. A typical picture of 
the pressure cycle is shown in Fig. 5. 

While the bubble chamber was warming up to its pre-expansion 
condition at 27.8°K and 70 psig, the Linlor gauge was calibrated by 
simultaneously recording the actual chamber pressure at 5-psig in- 
tervals as indicated on the Bourdon gauge, and photographing the voltage 
output on the scope for later measuring on a traveling microscope. 
Calibration data were also taken while the chamber was being cooled down 
at the end of the experiment. Data were also collected on other days in 
order to obtain some knowledge of the kind and amount of drift in the 
Linlor gauge. All the data, when analyzed, provided a reasonably good 
calibration curve of chamber pressure versus pulse height. 

The photographs of the Linlor-gauge output voltage taken during the 
experiment provided a measure of chamber pre-expansion and post- 


expansion pressure. With each photograph that was taken the actual 
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pre-expansion pressure was read from the Bourdon gauge, and com- 
parison of this pressure with the pulse height gave a measure of the 
Linlor-gauge drift for that particular chamber cycle. If the pre- 
expansion pressure and its corresponding pulse height did not fall on 
the calibration curve, the readings were normalized by use of the pre- 
viously obtained knowledge of drift. After normalization of the reading, 
the postexpansion pulse height as measured on the picture was applied 
to the calibration curve, and the corresponding pressure obtained. 

The value of the postexpansion pressure during the entire experi- 


ment was 48 + 5 psia. 
4. Thermodynamics of the Expansion 


As was mentioned in the section on temperature and pressure measure - 
ments, the pre-expansion and postexpansion pressures are known by 
‘direct measurement. The pre-expansion temperature is also known. 
The expansion of the liquid takes place in about 10 milliseconds, and 
therefore if there is a change in temperature during expansion itis not 
registered on the vapor-pressure meter, since the time constant of the 
meter and associated system is many times as large. There are two 
possibilities of importance; the expansion is either isothermal or 
adiabatic. The question of which process is correct was not answered 
experimentally, but has been calculated theoretically by solution of the 
boundary-value problem of heat transfer from the walls of the chamber 
to the liquid in the chamber [13]. The results show that the time for the 
temperature difference between the wall and the center of the chamber 
to drop to 1/e of its initial value is approximately 2 hours. The effects 
of convection were neglected in this calculation. The characteristic 
time involved in convection phenomena is certainly less than 10 milli- 
seconds, otherwise there would be large track distortions in the chamber. 
This provides rather convincing proof that the expansion is adiabatic. 

Knowing the expansion to be adiabatic, we can determine the post- 
expansion temperature of the liquid hydrogen. This, however, requires 
extrapolation from the liquid region into the superheated liquid region, 


using the latest appropriate experimental data on liquid hydrogen. 
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Johnston, Keller, and Friedman [14] found that-on a P-T diagram the 
isochores (lines of constant volume) are, within experimental error, : 
straight lines. Therefore extrapolation seems quite reasonable. 

The ratio of the adiabatic change of density with pressure to the 
isothermal change of density with pressure is equal to the ratio of CL 


t 


and c_., 1.e., 
p 


(9p/8P). € 


(9p/8P),,. 7c 


Vv 





p 


Extrapolating the data of Eucken [15] and Gutsche [16] as summarized in 
NBS Report No. 1932 [17] into the supe rheated region gives 3.2 cal/ 
mol-°K and 9.4 cal/mol- °K for the values of c,, and c_, respectively. 
This gives a ratio of 0.34. Using the pre-expansion values of T = 21 IK 
and P = 70 psig and the data of Johnston et al. [14], one finds the post- 
expansion temperature of the liquid naydrogen to be 27.5 + 0.1°K (a 
temperature drop of 0.3 + 0.1°K). | 

As was seen earlier, the pressure in the chamber is known to con- 
siderably less accuracy than the temperature. To see what the relative 
importance of these two quantities is in determining the density of the 


hydrogen one can find the ratio 
(8p/2T), fae) 
OT 
(9p/ 8P),. 


Johnson et al., [14] have represented their data with the empirical 


equation, 


P= Alp) - 7.01 T4437 pT , 


where FP is in atmospheres, p in sae) rte tT pin Ee, and A isa 


function of p only. Using this equation to determine (8F/ OT) at 
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T= 27.5°Kand P = 48 psia, one obtains 





(9p/ 8T) | 
ee = 57 = - 83,5 
(8p/ 8P),, K K 


Since our temperature measurements are accurate to 0.1 °K and our 
pressure measurements to 5 psi, we find the percentage error in density 
caused by the temperature error to be twice the percentage error caused 
by the pressure-measurement error. This then shows that the temperature 
is the thermodynamic variable that is more important in determining the 
accuracy of the density. | 

It was pointed out in the section on temperature measurements that 
there existed in the chamber a temperature gradient. With an average 
. vapor pressure difference of 3.5 psi (0.2°K), assuming the temperature 


gradient to be linear, we obtain a density gradient of 


dp /dz = 0.000059 g/cm*/cm 


and a percentage range gradient of 


(dR/R)/dz = 0.10%/cm. 


This gives a 1.6% change in range from top to bottom of the chamber. 
The effect of this on the average range of ur mesons in the chamber is 
to give a somewhat larger standard deviation than that caused by the 
3.3% expected range straggling. 

It should be pointed out that the data of Johnston et al. [14] are for 
liquid normal hydrogen (25% para-75% ortho). In the bubble chamber the 
liquid hydrogen differed in two respects. It was parahydrogen (99.70% 
para-0.21% ortho), and it contained a trace of deuterium. The content 
was as follows: 


H, = 99.03%, HD = 00.91%, D, = 00.06%, 
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Neglecting the D, content, the ratio of the range in parahydrogen to 
the range in the above mixture at the same conditions is dependent on 
the ratio of the specific volume of hydrogen to the specific volume of 


\ 


hydrogen deuteride in the following way: 


1 + (Waa) (V5. / V5.) 
ee 


ii R = eS tos a 
7) ad 1 + (%HD/%H,) 


Using the data from NBS Report 1932 [17], one finds this ratio to be 


Ru, / Rizp = 1.0009 + 0.0006 


Since this range error introduced by ignoring the deuteride content 
amounts to 0.1 + 0.1% of the final range, we considered it insignificant, 


and did not include it in arriving at our final results. 
5. System of Measurement of the ur Meson Track Length 


A. General 


As illustrated in Fig. 1, 35-mm stereo pair pictures were taken 
from above the bubble chamber. After it was developed, the film was 
scanned quickly for events of interest, in this case, a, ur sae 
decay chains. This initial scanning enabled us to eliminate those events 
which by their appearance in two dimensions obviously involved decay 
in flight of the a , or which by their orientation in the chamber could 
not be measured to the degree of accuracy desired. Those events which 
were not eliminated in this initial scanning were measured on a precision 
projection microscope developed at UCRL and shown in Figs. 6 and 7. 
The microscope projects the 35-mm pictures containing the events onto 
a viewing screen, The coordinates of any point in the pictures can be 
measured and punched automatically on IBM cards. This allowed us 
to store the coordinates of successive points along the tracks of 

+ 


7 an , abe e* particles on the cards, which when programmed into 


an IBM computer enabled us to determine the range of the vie” . The 


entire process is discussed in the following sections. 
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®. Brecision Projection Microscope 


The precision projection microscope shown in Figs. 6 and 7 projects 
the film pictures onto a large screen. The microscope table on which 
the film is placed is provided with independent x and y motions 
accurately determinable to the nearest micron. The position of the 
microscope table is transferred to IBM cards, using a Wang digitizer [18] 
and a standard IBM card-punch console. A fixed cross hair is also 
projected onto the screen and used in measuring points along the various 
particle tracks. The track length used for the t and e was approxi- 
mately one and a half times the length of the wp track. Seven to 12 points 
along each track were recorded on the cards in a definite sequence to 
simplify the computer program. The order in which tracks were 


measured is shown in Fig. 8. 
C. Solution of Optical Problem 


The space coordinates of a point in the chamber are obtained from 
the film coordinates by multiplying the film x and y coordinates 
with an appropriate magnification constant as follows (see Appendix I 
for details of the space construction): 

I I 
— = ee, y = kY, Zo = Wea. 
For a perfect optical system, the magnification constant k should be 
constant throughout any z plane, and be dependent on the horizontal 


lens separation (6) and the film coordinates of the point: 


em x 


The differences between the x, y, and z coordinates of two points in 
the chamber would then be 
I I 
x =Me, X ~ k., — ; 
= oa 
Ay k, Yy k, Y, ; 
AZ = (ky - k,) q . 
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Fig. 8. Schematic of stereo picture of bubble chamber showing 
the order of measurement of t> p> e events 
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The preceding remarks and those in Appendix I on the solution to the 
bubble chamber optical problem have been greatly simplified. Some of 
the complicating effects encountered and the methods used to correct or 
eliminate them should be mentioned. To determine k and 6, we had 
to measure on the film the distances between the images of reference 
fiducial marks permanently inscribed in the top and bottom glass. Then 
we compared these film distances to the known space distances between 
the fiducials. The results of these measurements required us to accept 
a different magnification for the x and y directions in any one z plane 
and a different magnification in each stereo view, or a total of four 
magnification constants for any one point. The need for these different 
magnification constants can be explained as the result of uneven shrink- 
age of the film during development, differences in the adjustments of 
the two lenses, anda preset 1.5° tilt of the chamber from the horizontal 
to prevent the accumulation of bubbles at the top glass. We added all 
the observed systematic variations in the optical system, and corrected 
for them by using different magnifications for each view and each 


coordinate direction within the view. 
D. Computer Program 


The problem of calculating the range of the Ta was programed for 
the IBM type 650 computer. Essentially the computer determined a 
least-squares parabola to fit each track in both views, and determined 
the points of intersection between the parabolas representing the 1 
and wt tracks and the a and e* tracks in each view. 

It is true that the track of a particle stopping in a magnetic field is 
not a parabola. But, weighing the ease of programing the problem for 
a parabola (as compared with a higher-order curve) against the required 
accuracy of measurement, we selected the parabola. Once the four points 
of intersection in the film were determined, it was a simple matter to 
have the computer determine the distance between the two intersections 
in space, as outlined in the section on the optical problem, as well as any 


other information desired about angles of intersection, average depth 
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of the pt iImethe chanieer, etc. To improve the over-all accuracy of 
the measurements, provision was made in the computer program to 
eliminate any event in which the coordinates of the last measured point 
at the end of a track was more than 9 microns off the least-squares 


parabola fit. 
E. Systematic Corrections to Measured Range 


We wish to determine the difference in length between the actual 
curved path (s) of the ut in the bubble chamber, and the straight-line - 
approximation (L) to the path calculated with the computer program. The 
approximate range-momentum relation 


3.66 (1) 


Ss = ap 


can be expressed in terms of the radius of curvature r in a magnetic 


field as 
oie eal 8 | (2) 


if the path of the ue is perpendicular to the magnetic field. To express 


this equation in x and y coordinates we use the parametric equations 


r= —; @im 7 =dy/de; ces 7 = dx/ds , (3) 


where 7 is the angle that a tangent to the particle path makes with the 


x axis. On substituting Eqs. (3) into (2) and using the relationship 


ae = suedere dr obtained by differentiating Eq. (2), we obtain 
re 3.668 , 2:66 | (4) 
2.66 
dye= sim eae suaeeue °° ae (5) 
ec leve: 
dx = Gq@s (aaa °2) sidaee dr. (6) 
2.66 | 
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Expanding the sin and cos terms in Eqs. (5) and (6) and performing the 
indicated multiplication, we obtain an infinite series. This series is 
convergent in the region of interest, that is, for particles of initial 
kinetic energy less than 5 Mev. Integrating the series term by term and 
summing significant terms, we can determine x and y accurately to 
four decimal places. We can then calculate L = (AS + ay?) if? and 

s-L 
solve for : 





, that is, the percentage difference between the length 
of the curved path and the straight-line approximation. 

We find this percentage difference to be 0.08% for the 1 and have 
ignored it in arriving at our final results. 

The mean range (R) as measured is not the true range (L) of the 
particle given by the Bethe-Bloch [3,4] theory. Because of multiple 
scattering, the measured range is shorter than the actual path length. 
To obtain an approximation of the difference between these two values, 
consider the particle path to be broken up into small segments dL, and 


the projection of these on the beam direction to be dR. Then we have 


dR = dLcos 0. , 


where 0. is the angle of scatter. Since we are considering small- 
angle scattering only, this gives 


q.¢ 
1 


dR = dL (1 - : 





Using the approximation by Segre [19] for the mean square angle of 


scatter, we have 


: dL 
age = dil -m,/m, In (E)/E.) = ; 


where Ep, is the initial energy and E, is the final energy of the y 


meson. With the approximate range-energy relation R= an ; 


we 
have integrated the above equation to obtain i: . Pieecorrec tien 


thus obtained is 0.1%. The value of range Sie in the results has been 
corrected by this amount. 
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6. Results and Conclusions 


During the experiment about 1400 stereo pair pictures were obtained. 
These yielded approximately 3000 1!’ ut +e? decay events in the 
chamber. During preliminary scanning we attempted to select for 
measurement only those events in which the angles of intersection of the 
particle tracks were between 30° and 150°. In addition, all events in 
which the ur track appeared to dip more than 45° were discarded be- 
cause of the large error inherent in the measurement of steeply dipping 
tracks. We were left with 771 events, which were measured on the 
precision projection microscope previously described. The computer 
program was designed to give not only the range of the ut but also its 
dip and the angles of intersection of the tracks. Of the 771 events 
measured, 327 hada ma dip angle greater than 30°, and were discarded 
as being subject to too large anerror. A histogram of the remaining 
444 events is shown in Fig. 9. Superimposed on the histogram is a 
normal distribution curve showing the average value of the range equal 
to the arithmetic mean of the 444 events and the standard deviation equal 
to that obtained from them. 

Taking the arithmetic mean as the most probable value, we obtained 
a range of 1.102 + 0.003 cm with a standard deviation of 0.055 cm (5.0 %) 
for the Te in liquid parahydrogen. The correction for multiple scattering 
gives 1.103 + 0.003 cm. This value corresponds toa z plane 7.1 cm 
from the top of the chamber, which is 16 cm deep. Using the range 
tables of Aron, [7], which are based on the Bethe -Bloch theory, we ob- 
tained the range-energy curves shown in Fig. 10. 

Using the compilation by Crowe [6] for rest mass and the initial 
energy of the ut meson (BE = 4.12 + 0.02 Mev and hae, = 206.86 + 0.l11m,), 
and the range tables of Aron[7], we obtained 0.0583 + 0.0006 igen for 
the density of superheated liquid parahydrogen at 27.6 + 0.1°K and 
48 +5 psia. The density under the same conditions for liquid normal 
hydrogen is calculated as 0.0591 ie” from Johnston's data [14], and 
0.0587 afl from Kammerlingh Onnes’ data [20]. It should be emphasized 


ae 





+ Abvoaileml Ge Spams 
ied pee ee LL ae a eee: ee 
od) 60@ ieee gee “ae “-"">» Ae | area iggt EAA vaekl 
Ve 10s od ees ee OOS Cee eg HOG «soup; 
SAE RR TAIN) Wee Wale aie Vine Heute Ghee Were eee 
oo eee UP aT co) 8 OF reeled eae pete! Ga 
aR treet ree OP) ced cot ee Geren dren % fl tone 
ring i retinte 4h mettre yr GF ca ey Scie gpl wih We, wes 
Or ee Fen Oierw dae tow IT) GUN ObSi eevee SW aden! 
kee AT ere eee ment CD ers Hole oe 
ee rememrreee™ GCN VHT WF GUNG HAL yl Gd Lule dow ner ews 
seme TTT wer ees ek oe te Rw EI fe baie Ue ine ih 
Sotrreerye? Gee Rage 0) cere seme oft vb "ye ued Tee pbevoceon: 
TNO RA ers a oyy) cee ae al Pulte geldc’ pe 
ee ee et ore Ls 
Oe a tee waite = lt Geta cere se (remeliseds Tearie 
Oe LOY He DCT WHOL Serer MAE GHEY HA SIO \f aa oy 
ce eS oe ne 
Cr err LOLS Ge Vad! 1+ concie 7 ree eit) fee 
OA oe PP) © ceed eee) Sarre BND Oe te) © oe ogee 6 
OV eke rel emer) EY osgorin etree diag om “oO srt pow 
we! ey ow a) TY me (Bod oe £41.) com 
eee a ee eee Geert 8) Re (ernee GAT pitt eft oven! 
SR et 8 it eet ection « FT) cmd be motes 
Tar were ep My e- wetwe 26d dead 
are 1 Kee feeo) mee tel fe) eed 99 charelip> sat geran 
te i ee 
re! ‘con DAA toe 8 80 8 eet er LT Long® fo REI wp et lous 
Cont a4 ote TP Le ree T eon apy. Te od (tins es 
Lemme Ew t) em rere thh Gih panty Mi aT tree # ome 
we ME) etd eee ee) ce gy leo pee af Qt 
el ee ee ee | fe a rn | 


e _> aa 
























*) 


RANGE HISTCGRAY yp MESONS 


62 Mean Ranga = 1.10220.002 em 
Standard Deviation = 0.055 em = 5% 


64 ( Number of » measured = 444 
60 
56 
52 


43 


49 
&] 36 


32 


NUMREP CF 


22 


20 


16 





TS | 


th 0SC*C:C“<CSCL CO Clk LIB Ioae 1036 130 Tae 
RANGE (em) 
MU-13398 


Fig. 9. Range histogram of 4.12-Mev a mesons in liquid hydrogen 
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Fig. 10. Range-energy curves inliquid H, bubble chamber based 
- on experimentally determined range of 1.103 4 .003cm for pt 
from mw decay. 
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that the latter two values are for normal hydrogen and that our value is 
for parahydrogen. Since the value of Johnston et al. [14] differs by only 
about 1%, we consider that their data extrapolated into the superheated 
region are adequate for bubble chamber thermodynamics and engineering 
purposes. 

Since there should be a vertical gradient in the chamber, we were 
interested in examining the correlation, if any existed between the 
range of the ur and its depth in the chamber. Using only those events 
(228) in which the ris track is very flat, 0° to 15° dip, we obtained a 
value of -0.14 for the correlation coefficient and determined the slope 
of the regression line of range on depth to be -0.002 cm/cm. This value 
compares favorably with that calculated by use of the temperature gradient 
observed and the data of Johnston et al., [14] which gave -0.001 cm/cm. 
The value obtained from correlation theory has a very large uncertainty, 
but does show that there is a variation of the range with vertical position 
in the chamber. It also points out that for particles that experience little 
straggling, the vertical position of the track should be taken into account 
when the energy is being calculated from residual range measurements. 

We have calculated the theoretical straggling by integrating 
numerically the expression given by Bohr [21], 
3 


0 
aR - R,)> ave. 4 eee, (dE/dx) ~ 


0 


= 3.2%, 


for the 4.12-Mev p* meson in liquid hydrogen. The value obtained is 
that amount of straggling caused by collisions with electrons only. To- 
ward the end of the fe track nuclear collisions would become important 
and give rise to a larger value. We have estimated our measurement 
error to be 2.5%, which with the 5% standard deviation obtained experi - 
mentally gives a straggling of 4.3%. We believe that this value represents 
a good estimate of the actual straggling of the 4.12-Mev a mesons in 


liquid parahydrogen. 
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APPENDIX I 


Solution to the Optical Problem 


Constants. Figure liis a diagram of the bubble chamber camera optics. 
No attempt was made to draw it to scale. For comparison purposes 
A= 100 cm, 6 =9 cm, q = 10cm. C is the center of the tap ghage of 
the chamber, D is the distance between the images of the chamber 
center as seen in the stereo pair. AF and A'F' are the optical axes. 


From geometrical considerations, we have 


AC/H=d/q «ad CA'/H’ = d'/q'; (A.1) 


assuming H =H', and q =q', we define the magnification in the plane of 
the top glass as k = H/q. Ce4), 
Substituting Eq. (A.2) in (A.1),we obtain 


AC =(d/q)H and) @A° = (d'/q)H; 
then 
AC +CA' =6& =(H/q) (d + d') 


§ 


k (D - 8) 


il 


6 ="D/tl = Way .. 
We can determine D and k directly from the film and then calculate 6. 
With these constants determined we can proceed to a determination of 
the coordinates of any point in the chamber from the stereo pair of 


pictures. 


Range Measurements in the Chamber. In the following refer to Fig. 1i. 


Point A is the assumed origin of coordinates in the bubble chamber, 
x,y, and z represent bubble chamber coordinates of point P, and 

a ae a, Tee are the corresponding film coordinates (subscripts 
denote point number, superscripts denote view number). The origin of 
the film coordinates is the projected optical axis in each view; positive 


X direction is to the right in the figure. We have the film coordinates 
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Fig. ll. Line drawing of 10-inch bubble chamber optical system 
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available from direct measurement; we must determine x, y, and z. 


Usine Mig. 11, then, we have 
J 
x/(H +z) = - X,/q 


and 
I 
(6 - x)/(H +z) =X, -/q 


solving these for H + z, we obtain 


eee (a/X,") =(8 - x) (q/+X)"), (A. 3) 


x =( [{H +2]/q) a eee ok, = 5/(X) - a), ee 


where ky is the mapnificationin tae z plane containing the point P. 
Similarily it can be shown that we have 


, I 
y=k,¥Y, 


Substituting Eq. (A.4) into (A.3), we get 


z=k,q-H. 


Finally, for the range of the o, we have 


Bao, __ ee 
Yr =" (es + AY +Az ) / : 
where 
A - & Gee - ee! 
14 PD 
I I 
Ay=k,¥, -k,Y, : 
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